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ABSTRACT 


Diffraction  radiation  is  that  electromagnetic  energy 
which  is  caused  by  a  relativistic  charged  particle  passing 
through  an  aperture  in  an  opaque  material.  Ter-Mikaelian 
solved  for  the  diffraction  radiation  from  a  point  charge. 
This  paper  discusses  the  phenomena  resulting  from  finite, 
relativistic  charge  bunches. 

Using  the  Huygens-Fresnel  principle,  diffraction 
patterns  from  spherical  and  cylindrical  charge  distributions 
are  found  and  plotted.  For  charge  bunch  sizes  less  than  the 
radiation  wavelength,  the  results  are  almost  identical  to 
those  for  point  charges. 

The  radiation  pattern  is  composed  of  two  regions.  The 
"transition  region"  is  characterized  by  a  strong  peak  at 
0  =  y”1,  the  Lorentz  factor.  The  "diffraction  region" 
consists  of  a  series  of  peaks  and  nulls  in  field  strength 
typical  of  the  standard  plane  wave  diffraction  pattern. 


4 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  -  7 

A.  BACKGROUND -  7 

B.  THE  VIRTUAL  PHOTON -  8 

C.  DIFFRACTION  THEORY  BASED  ON 

HUYGENS'  PRINCIPLE  -  10 

D.  FRESNEL  AND  FRAUNHOFER  DIFFRACTION  -  12 

II.  COMPUTATION  OF  DIFFRACTION  RADIATION  -  16 

A.  SOLUTION  FOR  A  POINT  CHARGE -  16 

B.  SOLUTION  FOR  A  SPHERICAL  CHARGE  BUNCH  -  18 

C.  DIFFRACTION  RADIATION  FROM  A  LARGE 

SPHERE  OF  CHARGE -  2  5 

D.  DIFFRACTION  RADIATION  FROM  A  CYLINDRICAL 

CHARGE  BUNCH  -  31 

E.  APPROXIMATIONS  AND  LIMITATIONS  -  38 

III.  PARAMETERS  AFFECTING  DIFFRACTION  RADIATION  -  39 

A.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  Y  -  40 

B.  DEPENDENCE  OF  DIFFRACTION  RADIATION 

ON  BUNCH  RADIUS -  4  3 

C.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON 

APERTURE  RADIUS  -  45 

IV.  DIFFRACTION  RADIATION  POWER  -  49 

V.  CONCLUSIONS -  51 

APPENDIX  A:  EVALUATION  OF  LINE  CHARGE  APPROXIMATION 

FOR  CYLINDRICAL  BUNCH  -  53 

APPENDIX  B:  DIFFERENCES  WITH  THE  THEORY  OF 

TER-MIKAELIAN  -  58 

LIST  OF  REFERENCES -  61 


5 


BIBLIOGRAPHY  - 

INITIAL  DISTRIBUTION  LIST 


6 


INTRODUCTION 


I . 


A.  BACKGROUND 

Diffraction  radiation  occurs  when  a  relativistic  charged 
particle  passes  through  a  hole  in  a  conducive  or  dielectric 
material.  This  phenomenon  is  very  well  understood  for  single, 
infinitesimal  charges.  Ter-Mikaelian  computed  the  diffraction 
radiation  for  a  single,  relativistic  electron.  [Ref.  1]  Less 
understood  are  the  effects  of  charge  bunching.  Presumably, 
bunch  sizes  on  the  order  of  the  radiation  wavelength  will 
significantly  affect  the  diffraction  radiation  intensity  by 
destructive  interference.  This  paper  discusses  diffraction 
radiation  from  charge  bunches  with  radii  greater  than  zero  but 
less  than  the  radiation  wavelength.  Diffraction  radiation 
will  be  computed  for  spherical  and  cylindrical  bunch  shapes. 

In  addition,  the  effects  of  bunch  size,  beam  energy,  and 
aperture  radius  will  be  explored. 

When  electrons  are  travelling  at  relativistic  velocities, 
they  can  be  treated  as  "pseudo-photons."  A  method  very  simi¬ 
lar  to  that  used  in  physical  optics  based  on  the  Huygens 
principle  is  then  used  to  find  the  diffraction  radiation  pat¬ 
tern.  [Ref.  l:p.  376]  Panofsky  and  Phillips  describe  the 
"virtual  photon  concept"  in  Reference  2. 
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B.  THE  VIRTUAL  PHOTON 


The  electric  field  from  a  "fast  electron"  approaches 
that  of  a  plane  wave  as  the  electron's  velocity  approaches 
c,  the  speed  of  light  in  a  vacuum.  Representations  of  the 
E-fields  of  electrons  at  various  velocities  are  shown  below: 


f 

m 


v  -  0  v  <  c  V  -  C 


For  a  uniformly  moving  electron,  one  can  solve  the 
Lienard-Wiechert  potentials  to  find  the  field  as  viewed  by 
a  stationary  observer.  One  finds  the  component  of  the  elec¬ 
tric  field  perpendicular  to  the  electron's  motion  is  very  much 
greater  than  the  parallel  component  and  is  given  by 


1 


eb 


4tts 


Y  ,.  2  ^  ,  .  >2,3/2 

o  (b  +  ( vty )  ) 


where : 


e 


the  electron  charge; 


b  =  the  perpendicular  distance  from  the 
electron's  path; 


Y 


/  2  -1 

the  Lorentz  factor  =  (  vl  -8  ) 


v/c . 
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One  can  view  this  field  as  a  sharp  pulse  of  electromagnetic 
energy.  The  magnitude  of  this  pulse  is  highly  dependent  on 
y,  which  defines  the  relativistic  energy  of  the  electron. 

In  addition,  the  duration  of  the  pulse  decreases  as  the 
electron  energy  rises.  The  electric  field  of  an  untra- 
relativistic  electron  will  appear  very  nearly  as  a  delta 
function  in  the  time  domain.  A  representation  of  this 
electromagnetic  pulse  is  shown  below: 


The  Fourier  transform  of  E^(t)  gives  the  frequency  components 
of  this  pulse. 


jgj1 


=  27  /  El(t) 


iut 


dt 


eby 

2 


-2 


iwt 


.  >  2  -2  2  3/2 

8 tt  eQ  -°°  [  (vt)  +y  b  ]  ' 


dt 


i  (— )  C 
vy  ^ 


8tt2£  bv  (£2  +1)  3/2 

o  ^ 


dC 
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where : 


r  =  IZZ 
4  b  ' 

Integrating,  Panofsky  and  Phillips  found 


Oil 


4  IT  ‘ 


£  bv 
o 


[  (— )K,  (— )  ] 
vy  1  vy 


where  is  the  modified  Bessel  function  of  the  second  kind. 

(x)  approaches  infinity  as  1/x  near  x  =  0  and  goes  to  zero 
as  e  '  for  arguments  greater  than  one.  Panofsky  and  Phillips, 
as  well  as  Ter-Mikaelian  in  his  development  of  diffraction 
radiation,  approximated  E  ^  by  equating 


E 


oo  1 


4tt 


e 

2  ebv 


for 


co 


< 


yv 

b 


0  for 


oo  < 


yv 

b  * 


In  his  treatment  of  diffraction  radiation  from  a  single 

particle,  Ter-Mikaelian  assumes  that  the  Fourier  components 

of  E  ^  are  spatially  limited  to  a  circle  of  radius  b  z  Ay . 

He  states  that  the  particle  will  radiate  at  a  frequency  oo 

only  when  Ay  >  r  ,  where  r  is  the  characteristic  dimension 

o  o 

of  the  aperture.  [Ref.  l:p.  377] 

C.  DIFFRACTION  THEORY  BASED  ON  HUYGENS'  PRINCIPLE 

Ter-Mikaelian  used  Huygens'  principle  to  compute  the 
diffraction  radiation  pattern  from  a  charged  particle  passing 
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through  an  opening.  He  thought  of  each  point  in  the  aperture 
as  a  source  of  secondary  waves  and  then  integrated  over  the 
aperture  to  find  the  total  field  at  some  distant  observation 
point.  The  field  at  each  point  within  the  aperture  is 
assumed  to  be  unaffected  by  the  obstruction.  In  practice, 
edge  effects  will  be  apparent  for  small  apertures. 

Huygens  Principle  states  that: 

every  point  on  a  primary  wavefront  serves  as  the  source 
of  spherical  secondary  wavelets  such  that  the  primary 
wavefront  at  some  later  time  is  in  the  envelope  of  these 
wavelets.  Moreover,  the  wavelets  advance  with  a  speed 
and  frequency  equal  to  that  of  the  primary  wave  at  each 
point  in  space.  [Ref.  3:pp.  60-61] 

Using  this  theory,  one  can  describe  the  propagation  of 

a  wave  past  an  obstruction  or  aperture  by  mathematically 

summing  the  contributions  of  each  "Huygens  wavelet."  Fresnel 

expanded  upon  Huygens'  principle  by  stating: 

every  unobstructed  point  of  a  wavefront,  at  a  given  instant 
in  time,  serves  as  a  source  of  spherical  secondary  wavelets 
(of  the  same  frequency  as  the  primary  wave) .  The  ampli¬ 
tude  of  the  optical  field  at  any  point  beyond  is  the 
superposition  of  all  of  these  wavelets  (considering  their 
amplitudes  and  relative  phases).  [Ref.  3:p.  330] 

There  are  some  problems  in  using  the  Huygens-Fresnel 
principle,  most  notably  the  effects  of  electron  oscillators 
at  the  edge  of  the  aperture.  Remember  that  Ter-Mikaelian 
assumed  that  the  field  within  the  aperture  was  unaffected  by 
the  obstruction.  Oscillating  electron  "clouds"  at  the  edge 
of  the  aperture  will  produce  their  own  contributions  to  the 
total  electric  field  within  the  aperture.  The  field  is  then 
made  up  of  a  part  from  the  incident  electric  wave  and  one  from 
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the  induced  oscillations  at  the  edge  of  the  obstruction. 
However,  for  apertures  larger  than  the  wavelength  of  the 
incident  radiation,  Huygens-Fresnel  diffraction  theory  is 
expected  to  provide  quite  accurate  results. 

D.  FRESNEL  AND  FRAUNHOFER  DIFFRACTION 

If  the  observation  plane  is  located  relatively  far  from 
the  aperture,  one  can  ignore  the  vector  aspects  of  wave 
propagation  and  use  the  scalar  Huygens-Fresnel  principles. 

For  each  small  source  at  the  primary  wavefront,  the  contribu¬ 
tion  to  the  field  at  observation  point  P(X,Y)  is  given  by 

dE  (X , Y)  =  E  (x,y)  exp  i (wt-kr)  ds/r 
P  ^ 

The  exponential  term  provides  the  phase  information  and 
the  r  term  in  the  denominator  takes  into  account  spherical 
spreading.  A  diagram  of  the  problem  is  shown  below: 
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In  the  typical  diffraction  problem  encountered  in  the 
study  of  physical  optics,  a  plane  wave  is  assumed  at  the 
aperture,  so  the  electric  field  across  the  aperture  has  no 
spatial  variation  in  amplitude,  polarization  or  phase.  When 
dealing  with  the  radiation  from  charged  particles,  however, 
one  must  consider  that  the  field  is  radial  with  amplitude 
changing  with  distance  from  the  charge.  In  our  problem,  we 
locate  the  observation  point  on  the  X-axis  and  deal  only 
with  the  x-component  of  the  field  at  the  aperture.  The  field 
is  symmetric  about  the  Z-axis. 

The  total  field  at  P(X,0)  is  given  by  the  superposition 
of  all  the  Huygens  wavelets  from  all  the  sources  at  P(x,y) 
of  size  ds  across  the  aperture. 

E  (X,0)  =  /  E  ( x , y )  e*P  i!f-kr)  ds 

Aperture 

This  is  a  very  difficult  integral  to  evaluate  in  that  the 
exponential  term  has  a  modulus  of  2ir ,  making  it  difficult  to 
make  an  asymptotic  approximation. 

The  phase  factor  can  be  approximated  by  taking  the  binomial 
expansion  of  r  and  approximating  it  to  the  first  order. 

9  o  9  999  22 

r  =  Z  +  (X-x)  +  (0-y)  =  Z  +  X  -  2Xx  +  x  +  y 

9  2  2  2 

=  (Z  +  X  )  -  2Xx  +  x  +  y 

=  R2 (1  —  2Xx/R2  +  (x2  +  y2 ) /R2 ) 
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where : 


2  2  2 
R  =  x  +  Z 


so 


r  =  R ( 1  -  2Xx/R2  +  (x2  +  y2)/R2)1/2 


The  binomial  expansion  is  given  by 


(1  +  x 


n  ,  ,  ,  n(n-l)  2 

=  1  +  nx  +  — — —  x  + 


To  first  order, 


r  =  R  -  Xx/R  +  (x2  +  y2)/2R 


The  total  field  at  P(X,0)  becomes 


E  (X,0I 
P/X 


=  e 


i  (wt-KR) 


Ea,x<x'y) 


Aperture 


x  ei^  (  (x2+y2)/2R  -  Xx/R)  ds 


For  large  R,  one  can  replace  the  r  in  the  denominator  with 
R  and  the  intearal  becomes 


E  (X,0) 
P'  X 


i  (ut— kr) 

5 - E -  /^X(x-y)e  R  ds 


14 


If  R  is  much  greater  than  the  aperture  size  and  the 

phase  difference  from  all  of  the  incremental  sources,  ds, 

2  2 

across  the  aperture  is  negligible,  then  the  k(x  +y  )/2R 
term  in  the  exponential  is  very  small  and  may  be  neglected. 
In  general,  if  this  term  is  less  than  one  radian,  then  the 
integral  can  be  approximated  by 

i  (wt-icR) 

E  (X,0)  =  - - = -  / E  (x , y ) exp ( ikXx/R) ds  . 

P  /  X 

In  the  Fraunhofer  approximation  the  electric  field  at  a 
distant  point  is  just  the  two-dimensional  Fourier  transform 
of  the  electric  field  in  the  aperture. 
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II . 


COMPUTATION  OF  DIFFRACTION  RADIATION 


A.  SOLUTION  FOR  A  POINT  CHARGE 

A  single  relativistic  electron  has  a  Fourier  spectrum, 
as  derived  earlier,  given  by 


ul 


4tt  ^eQbv 


i  <77>  Ki  <77' ) 


Using  the  Fraunhofer  method,  the  field  at  a  distant  observa¬ 
tion  point  is  given  by 


i  (to  t-«  R) 


E  (X,0,Z) 

P  ,  A 


Aperture 


_  ,  >  ikXx/R  JO 

Ea , x (x ' y ) e  dS 


where  the  x-component  of  the  field  at  the  aperture  is  given 
by 


Ea,  X 


4tt  ‘ 


£  bv 
o 


[  (— )K  A] 

Vy  1  Vy 


cos  0  . 


Changing  to  polar  coordinates  gives  us 


i  (cot-xR) 
e 

Z 


2  tt 

t  e  r  ,cob. ,mb\  1  _  ikXbcos9/R  ,  ,, 

— ~ -  ( — )K,  ( —  cos  6  e  b  d9db 

n  .  2  ,  vy  I  vy 

0  4tt  e  bv  1  ' 

o 


i  (cot-xR) 


K 


— ~  /2" 

4tt  £  v  0  0 


A  i(kbX/R)cos  9  ,, 

cos  0e  '  d0  db 
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Integrating  over  d0  yields 


which  contains  the  final  form  of  the  diffraction  radiation 
pattern.  As  Ter-Mikaelian  did  in  Reference  1,  the  radiation 
field  is  found  by  subtracting  the  particle  field  from  the 
results  above.  The  particle  field  is  found  by  letting  a 
go  to  infinity  in  the  integral  over  the  aperture.  It  is 
interesting  to  note  that  as  A  goes  to  infinity,  the  diffrac¬ 
tion  radiation  field  strength  goes  to  zero.  Physically,  at 
this  point  the  electron  no  longer  "feels"  the  presence  of  the 
obstruction . 

Ter-Mikaelian,  in  Reference  1,  computes  the  diffraction 
radiation  from  a  charged  particle  passing  through  a  circular 
aperture.  Using  the  method  outlined  in  the  previous  section, 
he  determines  (for  a A  <<  1) 


E 


p,X 


q  +  a 
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where : 


0) 

‘  ”  vy  ; 

q  =  k  sin  0  =  the  projection  of  the  wave  vector 

k  on  the  plane  z  =  0  ; 

K  =  the  aperture  radius. 

This  result  appears  to  be  significantly  different  from  the 
expression  derived  earlier.  In  reality,  however,  the  two 
plots  are  virtually  identical.  (See  Appendix  B.) 

B.  SOLUTION  FOR  A  SPHERICAL  CHARGE  BUNCH 

The  diffraction  problem  has  been  solved  for  single 
charged  particle  radiation. 

The  process  for  a  finite  charge  bunch  can  be  solved,  as 
a  first  approximation,  by  assuming  a  spherical  charge  dis¬ 
tribution  of  radius  aQ  with  constant  charge  density,  p. 

In  the  frame  of  the  moving  charge,  the  electric  field 
is  radial  with  a  magnitude  given  by 

for  r  <  a 

—  o 

for  r  >  a 

—  o 

This  result  is  easily  obtained  using  Gauss'  law.  A  repre¬ 
sentation  of  E(r)  is  shown  below. 
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To  find  the  field  from  the  moving  charge  distribution  in 
the  frame  of  a  stationary  observer,  one  must  find  the  Lorenz 
transformation  of  the  radial  field.  With  a  spherical  charge 
distribution,  an  observer  views  a  field  emanating  from  the 
center  of  the  charge  bunch;  the  field  from  the  uniform 
spherical  distribution  appears  to  be  a  point  charge  equal  to 
the  total  charge  of  the  sphere.  If  the  observer  is  located 
within  the  sphere,  he  still  "sees"  a  point  source,  but  of  a 
reduced  strength  (that  portion  of  the  sphere's  charge  enclosed 
by  a  Gaussian  surface).  Consequently,  the  Lorenz  transfor¬ 
mation  of  the  moving  spherical  charge  bunch  must  have  a  form 
similar  to  that  of  the  point  charge.  In  the  lab  frame,  one 
finds  a  perpendicular  electric  field  of 


(t) 


( 


4/3  tt  a  p 
o 


Y  2b 


4ttg 


4/3nb  p 
4tte 


r  .  .  ,  2  ^  -2.  2. 3/2 
[  (vt)  +  y  b  ]  ' 


Y  2b 


.  ,  _  2  ^  -2,  2,  3/2 
( (vt)  +  Y  b  ]  ' 


for  1  aD 


for  b  <  a 
—  o 
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This  form  of  E  (t)  is  only  approximate  for  b  <  aQ.  For 
small  bunches  (a  <<  XQ) ,  the  error  introduced  by  this 
approximation  is  negligible. 

Since  E  (t)  of  the  sphere  has  the  same  time  dependence 
as  that  of  the  point  charge,  it  must  have  Fourier  components 
of  the  same  form. 


E 


iw 


4/3Tra  p  ,  , 

_ _ oir  /Wb,  i 

.  2  .  vy  K1  vy 

4  tt  e  bv  1  ' 

o 


4/37Tb  p  r  (ub|  T,  ,wb,  n 
7~2~Zr_[  (vy)  K1  (vy)  J 


>  4tt  £Qbv 


-vy' 


for  b  >  a 
—  o 


for  b  <  a 
—  o 


Again,  E(  ^  is  negligible  for  A  <  b/y ,  we  make  the  additional 
requirement  that  the  wavelength  be  large  compared  to  the 
bunch  size  in  order  that  the  bunch  radiate  coherently. 

Again,  one  proceeds  using  the  Huygens-Fresnel  principle 
to  find  the  radiation  pattern  from  the  spherical  bunch  passing 
through  a  circular  aperture.  With  the  observation  point 
located  on  the  X-axis,  one  expects  to  see  only  an  X-component 
to  the  perpendicular  field.  Due  to  symmetry,  only  the  x- 
component  of  the  field  across  the  aperture  contributes  to 
the  total  field  at  P(X,0).  The  transverse  electric  field 
at  P(X,0)  is  given  by 


i  (lOt-KR) 


Ep/X(X,0,Z) 


Aperture 


„  ,  ,  ikXx/R 

Ea,x(x,y)e  dS 
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where  the  x-coraponent  of  the  field  at  the  aperture  is  given 
by 


E 


a ,  X 


4/37ia  p 
o 

4tt^£  bv 
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r  ,00  b  x  ,  (job  V  , 

[(— )K1(— )]c°s  e 
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for  b  <  a 
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for  b  <  a 
—  o 


Changing  to  polar  coordinates  yields 


E 


P/X 


i  (cot-KR) 


/Ea,X(b'6,e 


i  (kXbcos  0/R) 


ds 


where : 

ds  =  b  d Gdb  . 
Finally,  one  obtains 
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We  shall  first  solve  for  the  integral  from  b  =  0  to  a  ; 

o 


we  have 


/  °  f2”  VMA  [(ifc)K  (— )  ]  cos  0ei<te!tcos  0/E) 

0  0  4tt  ebv  ^  ^ 


b  d0db 


d6db  . 


The  integral  over  d0  is  equal  to  ( 2 tt i )  (kXb/R)  .  The  integral 
becomes 


2 tt i  f  b3(^)K,  (iafe)j,  (£^)  db 

4ti2  e  v  0 


1  vy 


R 


As  an  example,  one  can  expect  to  see  parameters  in  an 
experiment  on  diffraction  radiation  in  the  micro-wave  region 
as  follows: 


f 


a 

o 

y  . 
min 

v 


k 

X/R 


8  GHz 
0 . 5  cm 
20 

c  =  3  x  lO^  cm/s 
2tt/A  =  1.675  cm  ^ 

sin  0  where  0  is  the  angle  between  k  and  the 
Z-axis . 
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Using  these  values,  the  maximum  values  of  the  arguments  in 


(x)  and  J-^  (x) 

in  the  integral 

from 

0  to  a  are 
o 

(cob/yv) 

=  8.33  xlO-2 

(kXb/R) 

=  8.38  xlO-1 

With  these  values  for  arguments, 

the  integral  can  be  very  closely 

the  two  Bessel  functions 

approximated  by 

X 

1 — 1 

=  1/x 

and 


J1(x)  =  x/2  . 

The  integral  becomes 
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We  next  evaluate  the  integral  from  aQ  to  ftQ.  Since  the 
integral  here  has  the  same  6-dependence  as  in  the  previous 
case,  the  integral  becomes 
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Note  that  in  this  case,  we  have  not  made  the  approximations 
for  the  two  Bessel  functions  since  b  will  in  practice  extend 
to  values  much  larger  than  one  wavelength.  This  integral 
can  be  easily  solved  and  yields 


/  (“bjK-Aj.&db 

1  vy  1  vy  1  R 
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Combining  the  solutions  of  the  two  integrals  and  sub¬ 
tracting  the  particle  field,  the  total  transverse  field  at 
P(X,0)  is  given  by 


E  (X,0) 

P,  A 


3  2 
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C.  DIFFRACTION  RADIATION  FROM  A  LARGE  SPHERE  OF  CHARGE 

For  diffraction  problems  where  the  radiation  contribution 
from  within  the  charge  bunch  is  significant  (i.e.,  aQ  -  A  ) , 
one  can  not  make  the  approximation  used  in  the  previous  section 
where  E^  had  the  same  time  dependence  for  b  <  aQ  as  it  did  for 
b  >  aQ .  In  this  section,  we  shall  gain  a  more  exact  under¬ 
standing  of  the  problem  by  determining  the  diffraction  radia¬ 
tion  from  a  large  sphere  of  charge. 

Using  Gauss's  law,  one  finds 


( 


4/3,iao  b 


4  lie 


2  r 


for 


r  >  a 
—  o 


p  4/3Tir  b 
4tt£  2  r 


—  for 


r  <  a 
—  o 


where  b/r  provides  the  perpendicular  component  of  E,  as  shown 
below : 


V 
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This  relation  reduces  to 


E 

r 


P 

4  TT  £ 


4/3-rra^b 

3 

4/3  7T  b 

V 


for  r  >  a 

o 


for  r  <  a 
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A  representation  for  E^  is  shown  below  for  several  values  of 
b: 


This  function  for  b  <  aQ  can  be  integrated  numerically 
only.  However,  one  can  make  a  very  close  approximation  for 
E  with  b  <  aQ  by  writing 
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4/3Trb(1.5a  -Z) 

j  1.5aQ  -  'J a2-b2 

F1  =  4iF  ^  4/3"b 

4/3iTb(1.5a  +Z) 

1.5a  -  v42-b2 
o  o 
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o 

o 
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J  2  .2 
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o 

o 
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-  1.5a 

<  Z  <  -  n/ a2-b2 

o  o 


This  approximation  is  shown  on  Figure  1.  For  most  values 

3 

of  b,  this  linear  approximation  to  the  1/r  dependence  of 
E^ (r)  outside  the  sphere  lies  very  close  to  the  actual  curve. 
This  is  particularly  true  as  b  approaches  aQ ,  where  most  of 
the  diffraction  radiation  from  within  the  sphere  is  expected 
to  originate.  Replacing  z  with  -vt,  E  (t)  is  equal  to 


4/3TTaQb 
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for  b  >  a. 


12  2 

for  b  <  a  ,  v a  -b  <  -vt  <  1.5a 
o  o  o 


for  b  <  a  ,  -  '/a2-b2  <  -vt  <  '/a2-b“ 
o  o  o 


for  b  <  a  ,  -1.5a  <  -vt  <  - v a  -b 

o  o  o 


For  a  relativistic  electron  bunch,  the  Lorenz  transforma¬ 
tion  gives 
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The  diffraction  radiation  found  by  integrating  over  the 
aperture  from  aQ  to  fi  for  the  large  sphere  of  charge  will 
be  identical  to  that  found  in  the  small  sphere  approximation 
and  need  not  be  discussed  here. 
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This  section  will  discuss  only  that  radiation  originating 
from  the  aperture  for  0  <  b  <  aQ.  For  this  range  of  b,  the 
Fourier  components  of  the  field  are  found  by  writing  (t) 
as  a  difference  of  two  simple  functions  whose  Fourier  trans¬ 
forms  are  tabulated  in  any  introductory  text  on  Fourier 
analysis.  E^(t)  is  written  as  the  difference  of  two  tri¬ 
angular  pulses,  as  shown  below: 


It  is  easily  seen  that 
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A  diagram  of  a  typical  contribution  to  E  from  each  function 
is  shown  below: 


U 

One  can  easily  see  that  is  cut  off  at  high  frequencies 

The  field  at  a  distant  observation  site  P(X,0,Z)  is 
given  in  the  usual  manner  by 


E  (X, 0, Z) 
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lhis  integral  can  be  evaluated  numerically  to  find  the  dif¬ 
fraction  radiation  pattern  from  a  large  sphere  of  charge. 
When  this  is  done,  the  results  are  not  accurate  enough  to 
make  detailed  predictions  regarding  the  diffraction  pattern. 
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They  do  show,  however,  that  the  diffraction  pattern  from  a 
large  sphere  of  charge  is  consistent  with  that  found  for 
small  bunches.  The  field  strengths  are  significantly  less, 
but  the  pattern  of  peaks  and  nulls  is  roughly  the  same. 

D.  DIFFRACTION  RADIATION  FROM  A  CYLINDRICAL  CHARGE  BUNCH 
Consider  first  a  stationary  charge  distribution  with 
axes  as  shown  below: 


tC 


We  wish  to  find  the  electric  field  at  an  observation  point 
P (X , Y ) .  By  cylindrical  symmetry,  the  electric  field  at  an 

A  A 

observation  point  on  the  X-axis  will  have  only  i  and  k  com- 

/s 

ponents.  We  are  interested  primarily  in  the  i  component 
since,  for  the  relativistic  electron  bunch,  Ej^  >>  Ep  . 

In  order  to  solve  for  E^ ,  one  must  first  find  the  electro¬ 
static  solution  for  a  cylindrical  bunch.  One  can  make  an 
approximation  by  treating  the  cylindrical  bunch  as  a  line 
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For  b  <  a 


the 


charge  of  linear  charge  density  p(na^).  For  b  <  a  , 

2 

linear  charge  density  is  given  by  p  (Trb  ).  A  diagram  of  the 
problem  is  shown  below: 


For  b  >  a 
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For  a  slow  moving  bunch,  located  at  the  aperture  at  t  =  0, 
one  can  substitute  (-vt)  for  Z  and  write  the  time  varying 
field  at  the  aperture  as 
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For  a  relativistic  charge  bunch,  the  component  of  the 
electric  field  perpendicular  to  the  bunch  path  is  given  by 
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As  in  the  previous  section,  the  Fourier  components  of 
the  field  for  b  >  aQ  are  given  by 
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Evaluating  the  first  integral,  one  finds 
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letting  ip  =  (vty+D/b;  t  =  (b^/vy )  -  (L/vy )  ;  and  d^  =  (vy/b)dt, 
the  integral  becomes 
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Using  the  equation  e  =  cos  0  +  i  sin  0  and  remembering  the 
properties  of  even  and  odd  Fourier  integrals,  the  integral 
is  rewritten  as 
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Using  the  same  techniques  to  evaluate  the  second  integral 

in  E  ,  one  finds 
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As  previously  explained,  the  Hankel  transform  of  the  field 
at  the  circular  aperture  is  evaluated  to  find  the  field  at 
observation  point  P(X,0). 
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Simplifying,  the  equation  becomes 
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As  in  the  problem  for  a  spherical  charge  bunch,  integrating 
over  d9  gives 
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Evaluating  the  integral  from  0  to  aQ ,  one  makes  the  usual 
approximation 
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Next,  evaluate  the  integral  from  aQ  to  t  . 
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Since  X/R  =  sin  0,  the  total  field  at  the  observation  plane 
is 
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The  effect  of  bunch  length  is  in  the  (e  V^-l)  term.  Expand¬ 
ing  the  exponential  in  this  term  gives  an  imaginary  and 
real  part  of  the  solution. 
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At  first  glance,  this  function  seems  to  give  the  diffraction 
an  oscillating  dependence  on  radiation  frequency  for  all 
frequencies.  One  must  remember,  however,  that  the  K^(ojao/vy) 
term  in  the  solution  very  nearly  cuts  off  the  radiation  at 
a)  >  vy/aQ.  The  real  and  imaginary  terms  in  this  expression 
give  an  additional  phase  factor  to  the  overall  solution  for 
diffraction  radiation.  For  wavelengths  much  larger  than 
bunch  size,  the  real  term  is  approximately  zero  and  the 
dependence  on  bunch  length  is  given  by  i (— ) . 


E.  APPROXIMATIONS  AND  LIMITATIONS 

The  reader  is  cautioned  to  remember  that  in  making  the 
Fraunhofer  approximation  to  the  diffraction  problem,  one 
assumes : 

1)  The  diffraction  angle  is  small  (<  1  rad). 

2)  The  phase  difference  from  all  "sources,"  ds,  across 

the  aperture  is  small  (k(x^  +y2)/2R)  <  1). 

3)  The  aperture  size  is  greater  than  the  wavelength  of 
the  incident  radiation.  (Neglect  edge  effects.) 

4)  Since  charge  was  assumed  to  be  uniformly  distributed 
throughout  the  bunch,  the  quantized  nature  of  charge 
was  ignored.  This  can  be  done  if  the  electron  den¬ 
sity  is  sufficiently  high  (n-l/3  <  X)  . 

5)  Charge  bunch  sizes  are  assumed  to  be  less  than  the 
radiation  wavelength  to  ensure  coherent  radiation  of 
the  bunch.  Larger  bunches  have  smaller  intensities 
due  to  destructive  interference. 
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Ill . 


PARAMETERS  AFFECTING  DIFFRACTION  RADIATION 


In  the  previous  chapter,  diffraction  radiation  was 
computed  and,  for  a  given  frequency,  was  found  to  be  depen¬ 
dent  upon  several  parameters,  namely  y,  bunch  radius,  aQ , 
and  aperture  radius,  n.  .  In  addition  to  these  effects,  the 
radiation  is  significantly  cut-off  at  high  frequencies. 

Remembering  that  both  the  small  sphere  and  cylindrical 
charge  bunches  had  Fourier  components  proportional  to 
(aib/vy)  K(ojb/vy),  one  can  make  the  approximation  that  the 
fields  are  small  for  w  >  vy/b.  (When  the  argument  of  K^(x) 
is  1,  xK(x)  =  .6019.  This  gives  a  field  intensity  of  about 
one  third  the  maximum  value.)  If  the  maximum  value  of  b  is 
the  bunch  radius,  a  ,  then  the  maximum  value  of  oj  is  given 
by  ui  =  vy/a  .  One  expects  radiation  at  2oj  to  be  10  dB 
less  than  the  peak  values. 

Figure  2  can  be  used  as  a  rough  guide  to  the  frequency 
dependence  of  observed  radiation. 


X 


Figure  2 
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For  the  large  sphere,  the  Fourier  components  of  the 


fields  were  found  to  drop  off  as 


sin  (.75 


ooa 


vy 


U) 


One  can  say  that  the  field  is  cut-off  at  oj  =  vy/aQ. 

It  is  easily  seen  that,  for  each  bunch  shape  discussed, 
the  diffraction  radiation  is  very  near  the  peak  values  at 
microwave  frequencies.  A  previous  experiment  at  the  NPS 
Linac  discovered  diffraction  radiation  in  the  microwave 
region  (f  =  8  GHz)  from  a  circular  aperture  [Ref.  4].  For 
f  =  8  GHz,  ob/vyK^ (mb/vy)  is  very  nearly  equal  to  unity. 

The  plots  of  diffraction  radiation  vs.  0,  the  off-axis 
angle,  consist  of  two  distinct  regions.  The  first  is  charac¬ 
terized  by  a  strong  peak  at  9  =  y  \  Since  the  peak  at 
9  =  Y  is  typical  of  transition  radiation,  this  region  of 
the  curve  is  called  the  "transition  region." 

At  0  >>  y  ■*■ ,  the  radiation  pattern  appears  more  as  the 
typical  diffraction  pattern  of  peaks  and  nulls.  This  region 
is  called  the  "diffraction  region." 

A.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  y 

The  y-dependent  portion  of  the  function  describing  dif¬ 
fraction  radiation  tends  to  be  very  complicated  and  hard  to 
understand  analytically.  The  strongest  dependence  is  in  the 
transition  region,  where  a  very  strong  peak  in  field  strength 
occurs  at  0  =  y  ■*■ .  In  this  region,  the  field  strength  is 
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2  2 

very  nearly  given  by  (k  sinO)/((k  sin0)  +(to/vy)  )).  The 
strength  of  this  peak  is  directly  proportional  to  y.  In 
the  diffraction  region,  the  effect  is  less  pronounced. 

One  expects  to  see  more  diffraction  radiation  energy  from 
electrons  of  higher  energy.  Figure  4  shows  the  dependence 
of  diffraction  radiation  on  y.  The  curves  for  a  cylindrical 
bunch  will  be  identical  to  this  figure.  Increasing  y  tends 
to  sharpen  the  diffraction  pattern,  increasing  the  strength 
of  the  peaks  in  field  strength. 

For  pseudo-photons  of  a  given  frequency,  the  transverse 
fields  tend  to  extend  farther  out  for  increasing  y.  Alterna¬ 
tively,  the  field  for  a  given  frequency  and  radial  distance, 
b,  tends  to  be  higher  for  an  increasing  y.  Figure  3  demon¬ 
strates  this  point.  If  one  compares  the  same  Fourier  com¬ 
ponents  of  two  relativistic  electrons  with  different  energies, 
the  more  energetic  electron's  fields  extend  much  further  out 
than  those  of  the  less  energetic  electron. 


Figure  3 
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ripure  4:  E  vs  0  for  various  enerpies 


Ter-Mikaelian  determined  that,  for  a  point  particle 
with  cotQ/vy  <<  1,  the  y-dependence  was  given  by 

g  _ k  sin  9 _ 

(k  sin  0)2  +  (— )  ^ 
vy 

The  y-dependence  found  in  the  last  chapter  for  both  cylin¬ 
drical  and  spherical  charge  bunches  is  consistent  with  that 
of  Ter-Mikaelian.  One  expects  the  fields  to  go  to  zero  for 
non-relativistic  particles. 

B.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  BUNCH  RADIUS 

Throughout  this  paper,  we  have  assumed  that  the  bunch 
radius  is  smaller  than  the  wavelength  of  the  radiation.  When 
this  is  the  case,  one  may  assume  coherent  radiation  from 
the  bunch. 

By  examining  the  equations  for  diffraction  radiation  for 
all  three  bunch  types,  one  can  see  that  diffraction  radiation 
tends  to  decrease  for  bunches  with  larger  radii.  For  small 
bunch  radii,  a  general  decrease  in  the  field  strengths 
from  the  peak  levels  for  infinitesimal  bunch  sizes  is  observed, 
as  shown  in  Figure  5.  For  cylindrical  bunches  of  equal 
charge,  the  diffraction  radiation  field  strength  is  inversely 
proportional  to  L,  where  L  is  the  bunch  length.  This  is 
expected  since  the  radial  fields  at  the  aperture  are  propor¬ 
tional  to  the  charge/length  within  the  bunch.  The  (ewL//v^-l) 
factor  also  specifies  a  length  at  which  coherent  radiation 
from  the  bunch  can  be  expected. 
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0  (RfiDlftV s) 

Figure  5:  E  vs  9  for  various  bunch  radii 
3f=100,  R  =  10  cm  ,  f=  8GHz 


One  can  easily  understand  the  effect  of  bunch  size  by 
examining  the  static  radial  fields  from  a  sphere  of  charge. 
The  effect  from  a  cylindrical  charge  bunch  will  be  similar. 


Figure  6  compares  the  fields  from  two  spheres  of  radii 

aQ  and  a^.  As  shown  earlier,  for  an  infinitesimal  sphere, 

the  fields  increase  assymptotically  to  infinity  as  r  goes 

to  zero.  For  a  sphere  of  radius  a  ,  the  field  increases 

o 

with  the  same  1/r  dependence  as  r  decreases  from  infinity 
to  ac  and  then  drops  off  linearly  to  0  at  r  =  0 .  One  can 
easily  see  from  the  figure  that  the  sphere  of  smaller  radius 
is  the  source  of  more  electrostatic  energy  than  the  larger 
sphere.  As  the  sphere  gets  smaller,  the  total  electrostatic 
energy  rapidly  approaches  that  of  a  point  charge. 

C.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  APERTURE  RADIUS 
A  plane  wave  incident  upon  a  circular  aperture  results 
in  the  familiar  Airy  diffraction  pattern.  The  large  central 
maximum  in  intensity  is  known  as  the  Airy  disk.  The 
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diffraction  radiation  pattern  from  a  relativistic  electron 
beam  is  very  similar,  except  that  the  central  maximum  is 
replaced  by  a  null.  This  central  null  results  from  the 
radial  nature  of  the  electron  fields.  In  the  Fraunhofer 
diffraction  method,  integrating  over  the  aperture  for  d0 
results  in  a  Oth  order  Bessel  function  of  the  first  kind 
for  the  plane  wave  and  a  1st  order  Bessel  function  for  the 
electron  diffraction  radiation,  thus  explaining  the  princi¬ 
pal  difference  in  the  two  radiation  patterns. 

Figure  7  shows  the  diffraction  radiation  pattern  for 
cylindrical  and  small  spherical  charge  bunches  passing 
through  circular  apertures  of  various  radii.  Changing  the 
aperture  size  affects  the  diffraction  radiation  pattern  in 
both  the  transition  and  diffraction  regions.  One  can  define 

-+  (jj  ^ 

a  dimensionless  parameter,  k A  ,  where  k  =  — k  =  the  wave 
vector  and  A  is  the  aperture  radius.  This  parameter  alone 
determines  the  location  of  the  peaks  in  the  diffraction  region. 
When  ktQ  increases,  the  peaks  gain  in  strength  and  move  in 
closer  to  the  0  =  0  axis.  Decreasing  k-t  spreads  the  dif¬ 
fraction  pattern  out  over  wider  angles.  This  is  predicted 
since  the  diffraction  radiation  pattern  is  caused  principally 
by  the  ^  (k-tQsin  0)  terms  in  the  solution.  Note  that 
patterns  with  a  large  number  of  maxima  and  minima  are  found 
only  for  aperture  sizes  much  larger  than  the  wavelength  of 
the  observed  radiation.  This  is  to  be  expected  since,  for 
the  normal  diffraction  of  a  plane  wave,  a  small  aperture  acts 
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Figure  7:  E  vs  9  for  various  aperture  radii 
a  =  .1cm,  100,  f  =  8  OHz 


almost  as  an  isotropic  point  source  of  radiation.  In  the 
case  of  diffraction  radiation,  we  still  see  the  null  at 
8=0  due  to  the  radial  nature  of  the  fields.  Since  the 
dimensionless  parameter  is  given  by  k-t  ,  one  can  change  the 
appearance  of  the  diffraction  pattern  by  changing  fi  ,  the 
physical  size  of  the  aperture,  of  k,  which  determines  the 
wavelength  of  the  radiation. 

In  the  transition  region,  changing  the  parameter,  k ti  , 
affects  the  strength  of  the  peak  at  0  =  y  \  As  k/i  goes 
to  zero  the  peak  at  y  ^  increases,  asymptotically  approaching 
the  value  characteristic  of  transition  radiation.  As  k x 

o 

increases,  the  strength  of  this  peak  decreases,  until,  as 
k4Q  goes  to  infinity,  the  peak  goes  to  zero.  At  this  point, 
the  charge  no  longer  senses  the  obstruction  and  no  energy  is 
radiated . 

One  should  also  notice  that  the  diffraction  radiation 
pattern  for  spherical  and  cylindrical  charge  bunches  are 
identical  for  bunch  lengths  less  than  the  radiation  wave¬ 
length.  For  very  large  charge  bunches,  interference  effects 
are  expected  to  reduce  the  overall  intensity  of  radiation/ 
(bunch  charge) ,  but  the  location  of  maxima  and  minima  will 
not  change.  The  locations  of  the  radiation  pattern's  maxima 
and  minima  are  principally  a  function  of  radiation  wave¬ 
length  and  aperture  size,  as  it  is  in  the  normal  plane  wave 
diffraction  problem. 
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IV.  DIFFRACTION  RADIATION  POWER 


It  has  been  suggested  that  diffraction  radiation  may 
be  used  for  electron  beam  targetting  or  in  a  non-interceptive 
beam  diagnostic  device.  [Ref.  5] 

Present  diagnostic  systems  are  based  on  the  beam  passing 
through  a  thin  foil  and  using  transition  or  Cerenkov  radia¬ 
tion  to  analyze  beam  characteristics.  The  problem  with  this 
scheme  is  that  the  beam  is  disturbed  by  the  foil.  If  the 
beam  passes  through  an  aperture,  one  can,  in  theory,  obtain 
the  same  diagnostic  information  without  disturbing  the  beam. 
If  diffraction  radiation  is  to  be  used  as  a  diagnostic  or 
targetting  means,  one  must  gain  an  appreciation  for  the  power 
levels  expected  from  a  typical  device. 

For  an  electromagnetic  wave  given  by 

E  -  Eoei(“t-KR>  , 

the  power  per  unit  area,  or  irradiance,  is  given  by 

CE  - 

I  5  <S>  E  E  . 

2  o 

One  can  compute  the  diffraction  radiation  irradiance  by 
using  this  expression  together  with  the  equations  derived 
in  Chapter  II. 
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As  an  example,  the  diffraction  radiation  power 
delivered  by  a  1000  amp,  50  MeV  electron  beam  passing 
through  a  20  cm  diameter  circular  aperture  will  be  maximized 
at  0  =  10'.  If  one  assumes  10  cm  long  cylindrical  bunches 
of  radius  .1  cm  and  evaluates  the  expression  for  diffraction 
radiation  from  a  cylindrical  charge  bunch,  the  irradiance 
at  Z  =  1  m  from  the  aperture  is  equal  to 

CE  ~  9 

I  =  =  500  KW/m  . 

The  power  per  unit  solid  angle  is  500  KW/SR. 


50 


V.  CONCLUSIONS 


In  this  paper,  the  equations  regarding  diffraction 
radiation  from  a  relativistic  charge  distribution  passing 
through  an  aperture  have  been  derived  from  the  Huygens- 
Fraunhofer  treatment  of  the  fields  across  the  aperture.  The 
calculations  are  very  similar  to  the  familiar  problem  of 
plane  wave  diffraction  solved  in  every  physical  optics 
class,  with  the  complications  that  the  fields  are  radial 
with  varying  field  strength  in  space. 

While  a  plane  wave  incident  on  a  circular  aperture  causes 
a  diffraction  pattern  with  a  large,  central  peak  known  as 
the  Airy  disk,  the  diffraction  radiation  pattern  produced 
by  electrons  is  characterized  by  a  null  at  0  =  0  and  is 
highly  dependent  upon  aperture  size.  In  addition,  a  strong 
peak  at  0  =  y  is  present,  characteristic  of  transition 
radiation.  If  the  aperture  size  is  very  small,  the  structure 
in  the  diffraction  region  disappears  and  the  radiation 
pattern  approaches  that  of  pure  transition  radiation.  The 
only  radiation  present  is  the  strong  peak  at  0  =  y  .  In¬ 
creasing  the  aperture  size,  or  more  properly,  the  dimensionless 
parameter  k h.  ,  decreases  the  peak  at  0  =  y  ^  and  compresses 
the  structure  in  the  diffraction  region.  Eventually,  as  kA 
approaches  infinity,  the  radiation  goes  to  zero.  The  diffrac¬ 
tion  radiation  pattern  in  the  diffraction  region  exhibits 
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some  dependence  on  electron  energy  but  not  nearly  as  strong 
as  in  the  transition  region.  The  radiation  intensity  in 
both  the  transition  and  diffraction  regions  has  a  minor 
dependence  on  bunch  size. 

Ter-Mikaelian  describes  diffraction  radiation  as  having 
a  strong  peak  at  0  =  y  ^  followed  by  a  broad  series  of  peaks 
and  nulls  dependent  on  Jq  (k  sin04  ).  The  diffraction  pat¬ 
terns  developed  in  this  paper  are  primarily  dependent  on 
(k  sin0^o)  and  (k  sin  @4q),  but  when  plotted  are  found 
to  be  in  complete  agreement  with  Ter-Mikaelian ' s  results. 
Both  results  are  in  agreement  on  several  key  aspects  of 
diffraction  radiation: 

1)  There  is  a  strong  dependence  of  diffraction  radiation 
on  aperture  size. 

2)  The  diffraction  radiation  energy  goes  to  zero  as 
aperture  size  goes  to  infinity. 

3)  The  diffraction  radiation  energy  is  dependent  on  the 
energy  of  the  beam,  described  by  y. 

4)  The  diffraction  radiation  exhibits  a  strong  peak  at 

0  =  y"l. 
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APPENDIX  A 


EVALUATION  OF  LINE  CHARGE  APPROXIMATION 
FOR  CYLINDRICAL  BUNCH 

As  described  in  Chapter  II,  the  static  fields  for  a 
cylindrical  charge  distribution  were  approximated  by  a 
modified  line  charge.  The  linear  charge  distribution  was 
given  by 


A 


pTiao  C/cm  , 


where : 


a 

o 

P 


cylindrical  bunch 

charge  density  of 
bunch  in  C/cm-3 . 


radius  in  cm 
cylindrical  charge 


When  b,  the  off-axis  distance,  is  less  than  aQ  the  linear 

2 

charge  density  was  given  by  A  =  pub  . 


X 


For  the  geometry  shown  above,  the  field  for  the  line 
charge  is  given  by 
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p  7T  a 

i  (. 

4tt£  b  ' 


Z-L 


(b2  +z2)1/2 


(b2  + (Z-L) 2) 1/2 


As  a  check  of  this  approximation,  the  fields  predicted 
by  the  linear  approximation  were  compared  to  the  fields 
resulting  from  the  actual  cylinder  of  charge,  which  were 
computed  by  numerically  integrating 


P 

4tte 


/ 

r=0 


a 

o 


_ X  -  r  cos  9 _ 

( (X-r  cos  0 )  2+  (r  sin  0 ) 2  +  (Z-Z')2)3//2 


r  drd0  dZ 


Figure  A. 1.1  shows  the  comparison  for  three  "long" 
cylinders  (L  >>  aQ) ,  where  the  observation  point  is  outside 
of  the  cylinder  (b  >  aQ) .  The  linear  approximation  is 
expected  to  work  best  in  this  situation. 

As  can  be  seen  from  the  figure,  the  approximation  is 
very  close  for  the  cylinders  of  lengths  1  cm  and  2  cm  and  is 
identical  for  the  longer  10  cm  bunch.  The  radial  fields 
become  very  nearly  equal  to  that  of  an  infinitely  long  line 
of  charge  very  close  to  the  ends  of  the  bunch. 

The  linear  approximation  is  expected  to  perform  worst 
for  a  relatively  short,  thick  cylinder  of  charge  when  the 
observation  point  is  on  the  interior  of  the  cylinder. 

Figure  A1 . 2  shows  the  comparison  for  these  cylinders 
ranging  in  relative  length  from  L  =  aQ  to  L  >>  aQ.  For  the 
shorter  length,  the  approximate  fields  are  found  to  rise 
slightly  more  sharply  and  to  higher  levels  than  the  actual 
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Figure  8:  E  (‘ti7t/n  )  vs  Z  for  cylinders  of  various  lengths 
art =  1  cm,  b=  .5  cm 
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fields.  Even  in  this  most  extreme  example,  however,  the 
approximation  should  be  adequate  for  our  purposes.  As  the 
bunch  length  increases,  the  actual  and  approximate  fields 
are  found  to  agree  very  closely,  with  a  slightly  steeper 
rise  in  the  linear  approximation. 

In  conclusion,  we  find  that  treating  the  cylindrical 
charge  distribution  as  a  modified  line  charge  offers  a  very 
satisfactory  approximation  to  the  actual  fields  present. 

Using  this  approximation  for  a  static  charge  bunch,  one  can 
compute  very  closely  the  Fourier  components  of  a  relativistic 
charge  bunch. 
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APPENDIX  B 


DIFFERENCES  WITH  THE  THEORY  OF  TER-MIKAELIAN 


In  Chapter  III,  it  was  brought  to  the  reader's  attention 
that  the  equations  derived  in  this  paper  were  slightly  at 
odds  with  the  results  presented  by  Ter-Mikaelian  in  Reference 
1.  In  his  paper,  Ter-Mikaelian  stated  his  results  for  an 
electron  passing  through  the  center  of  a  circular  aperture 
of  radius  aQ  for  wao/vy  <<  1  (in  Gaussian  units) . 
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JQ  (qtQ)  cos  ip 


Jo  (qtG)  sin  ip 


where : 


a  =  w/vy 

q  =  k  sin  8  =  projection  of  the  wave  vector  on 

the  z  =  0  plane 

ip  =  the  angle  between  q  and  the  X-axis 


These  equations  give  a  radial  field  equal  to 
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J±  n  2  2 

2tt  c  q  +  y 


2  Jo(qV 


Since  this  equation  appears  to  be  inconsistent  with  the 
equations  derived  in  this  paper,  this  appendix  will  explain 
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at  length  the  method  used  to  obtain  the  Bessel  function  in 
the  equation. 

The  integral  over  d0  in  the  Fraunhofer  diffraction  problem 
was  written  as 

2  7T 

(  .  i  (kbX/R)  cos  0  ,Q  ,, 

}  cos  0  e  '  d6  db  . 

0 


The  cos  0  term  in  this  integral  is  due  to  the  radial  nature 

1  6 

of  the  electric  field.  Using  the  expression  cos  0  =  (e  + 

_  T  Q 

e  ) /2 ,  the  integral  can  be  re-written  as 


,e^  +  e  i  (kbX/R)  cos  0 

/  ( - 2 - >e 


d0 


-  T  / 


2 tt  i  (0+^^-cos0  )  i  (-9+ -^-cos8  ) 
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+  e 


d0 


T  ,  .  l  r  mv+ucosv  , 

Since  J  (u)  =  — ^ —  e  dv,  or 

m  2  tt  o ; 


2  jf 

~  .m_  ,  .  f  i (mv+ucosv)  dv  , 

2tti  J  (u)  =  e 

m  o 


the  integral  over  d0  is  equivalent  to 


TTiJ.  (kXb/R)  +  tt i  (kXb/R) 

1  -1 


.  ,  _  ,kXb<  T  ,kXb,  . 
lrl  J1  ~ R~  ‘  J-l(— >> 
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Since  J'(x)  =  {j  ,  (x)  -  J  ,(x)}/2,  with  n  =  0,  this 
n  n-i  n+i 

expression  reduces  to 


-2iTi 


Jo( 


kXb 

R 


which  equals 


o  •  t  A^b, 
27,1  Ji(— 1 1 


since  CP(x)  =  -J-^  (x)  . 

As  shown  earlier,  the  sin  8  dependence  of  the  diffraction 
pattern  is  given  by 

2  ~ — j f^0^Ki  )  J2  _  aJl  (q^0)  k2  <'a'x q)  )  J  • 

q  +  a 

This  appears  to  be  quite  different  than  that  found  by  Ter- 
Mikaelian . 

However,  if  one  plots  the  two  solutions,  he  will  find 
they  are  identical. 
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